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SUMMARY 


Various features of the solutions of Duffing’s equation are described using a representation of the 
solutions in the Laplace-Borel transform domain. An application of this technique is illustrated for the 
symmetry-breaking bifurcation of a hard spring. 


INTRODUCTION 


For linear systems, it is well recognized that the operational calculus based on Laplace transfor- 
mation is a convenient and powerful technique. One would like to have a similar operational calculus for 
treating nonlinear systems. Fliess and his coworkers (ref. 1) (cf. also ref. 2) advocate the use of functional 
expansions based on noncommutative power series. One practical consequence of the approach is the use 
of Laplace-Borel transformation in a similar way to that of ordinary Laplace transformation. In addition 
to the advantages of the Laplace transformation, the new transformation enables treatment of the product 
operation of functions via the so-called shuffle operation or melange. 

In this paper, we shall apply Laplace-Borel transformation to Duffing’s equation. This simple 
system shows some rich features of nonlinear mechanics, e.g., existence of symmetric and asymmetric 
regimes, transitions into chaotic regimes via an infinite number of period doublings. 
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The main difficulty in analytic approaches is the representation of solutions. Currently, we observe 
that most researchers are using iterative procedures and harmonic balance techniques. We shall discuss the 
validity of these approaches by exploiting the symmetry properties of Duffing’s equation. By applying the 
Laplace-Borel transformation, we shall show that the first-order approximation of the iterative solution will 
determine the nature of the solution. Given the first-order approximation, one is still faced with algebraic 
manipulation of nonlinear terms. Traditional harmonic balancing results in unwieldy calculations when we 
want to take into account the higher-order approximations. The Laplace-Borel transformation approach 
improves the harmonic balancing method with a more manageable algebra. We shall show that the Laplace- 
Borel transformation approach provides an operational calculus similar to that of the Heaviside calculus. 

Next, we shall apply this new algebraic technique to elucidate an intricate feature of the nonlin- 
ear mechanics of a hard spring, namely, a symmetry-breaking bifurcation. In this case, we shall show 
that only a second-order harmonic balah£6 approximation is sufficient to capture the essential features 
of a symmetry-breaking bifurcation. The bifurcation analysis based on Floquet’s theory gives values of 
parameters in reasonable agreement with exact numerical results. 


ALGEBRA OF LAPLACE-BOREL TRANSFORMATION 


We shall give a brief summary of the principal rules of Laplace-Borel transformation and refer the 
interested reader to references 1-3 for greater detail. To each analytical function f(t) having an expansion 

00 t n 

/<»- E%r (1> 

w=0 u - 

there exists a corresponding generating power series G given by: 

S=fwi (2) 

n=0 


We have the following algebraic rules which are obtained in references 1-3: 



<=> Xq 


(3) 


The Laplace-Borel transformation operates as if the integral is replaced by the symbolic variable x 0 ; while 
in the case of the Laplace transformation, the derivative is replaced by the symbolic variable s. The 
Laplace-Borel transformation of elementary functions is analogous to the Laplace transformation and we 
shall use the following relation: 

-jwxo)- 1 (4) 

The Laplace-Borel transformation has the additional capability to handle the product of functions via the 
so-called shuffle product 13 

f(t).g(t)^G(f)DG(g) (5) 
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where 13 denotes the shuffle product operation. For instance, the shuffle product of two rational simple 
functions is: 

( 1 - ox 0 ) _1 0 ( 1 - bx o)" 1 = ( 1 - (o + b)x 0 )~ l (6) 

The Laplace-Borel transformation treats the forcing term as a symbolic variable x\ 

[ f(r)dT<=>x i (7) 

Jo 

where / denotes the forcing term. The product of xi by xo is noncommutative 


xi.xo^xo.xi ( 8 ) 

because: 

f dri f f(n)dnj f f(7\)dri f dri (9) 

Jo Jo Jo Jo 

At the end of symbolic operations, x\ is replaced formally by: 

xi.=>x 0 .G(f)U (10) 


SYMMETRY PROPERTY OF DUFFING’S EQUATION 


Duffing’s equation is given as 

x + ax + wlx + fix 3 = fcoswt (11) 

where a > 0 and 0^0 . Duffing’s equation has an invariance property with respect to the following 
transformation: 

it 

t => t + — x =► — x (12) 

w 

Due to this property, a symmetric solution x, is defined as: 


*,(*) = -*,(*+-) (13) 

iv 

If a symmmetric solution is represented by a Fourier series, the series has only odd terms due to the property 
of equation (13), i.e., 


oo oo 

x,(t) = 52A2’*icos((2n+ + !)wt) (14) 

n=0 n^O 

The symmetric solution has a zero time average. By contrast, the asymmetric solution will contain even 
terms and hence may have a non-zero time average. 
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ANALYTICAL REPRESENTATION OF SOLUTIONS 


Let us integrate Duffing’s equation twice: 

x(t) + a Jo x(ri)dri + t v% fo dn Jo x(n)<fo + Pfo dn J 0 n x*(n)dn 

= f dn f dnf(n) + x( 0) + (x(0) + ax(0)) f dry (15) 

Jo Jo Jo 

where f(t) = fcos( wt ) and x(0) , x(0) are initial displacement and velocity respectively. Next, we shall 
obtain the Laplace-Boiel transformation of this equation by means of the algebraic rules listed earlier: 


G(x) + atxoGix) + wlxlG(x) + Px 2 0 G(x)l3G(x)UG(x) 

= xoxi + xo + (x(0) + ax(0))xo 
After some algebraic manipulation, we have 


(16) 


G(x) = -p(l -oix 0 )- 1 x 0 (l - a2Xo)- x x 0 G{x)UG{x)UG{x) 

+ ( 1 — OiXo) -1 Xo(l — 02 Xo) -1 Xl + (1 — OiXo) -1 ( 1 — 02Xo) -1 (x(0) + (x(0) + ax(0))Xo) (17) 

where af 1 and oj 1 are roots of the equation: 

1 + axo + tujjXo = 0 ( 18) 

Our interest lies principally in the large-time behavior of the solution. If we neglect the transient terms 
(i.e., terms that die off with large time) we have to the first order 


(^(xo) * (1 — oixo) -, xo(l — <X2Xo)~ x xoG(fcoswt) 
or, using the algebraic rule (eq. (4)), 

G ,(1) (x 0 ) = (1 -oix 0 ) _1 x 0 (l - 02Xo) -1 xo(//2)(l - jwx o) -1 + (1 + ;tux 0 ) -1 ) 
If we define 


0_i = 


1 


Wq —w 2 + jaw 
1 

*Ll “ 2 2 

Wq — W i — jaw 

then after decomposition into simple rational fractions and neglecting the transient terms: 
G (,) (x o) * (//2)(fl_i(l —jwx 0 )~ l + 0 , 1(1 + ;tux 0 )“ 1 ) 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 
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If we insert G (1) (xo) into the shuffle product term, we get a second-order approximation. The shuffle 
product is an operation which can be done on the computer using any one of the symbolic programming 
languages available to the analyst. In this analysis, however, we shall simply use the algebraic rule (eq.(6)) 
to obtain: 


G13G13G = 3Q^!f2+i( 1 — jwxo ) -1 + 3Q_i£2+,( 1 + jwx o ) -1 

+ Qi,(l — 3jwxo)~ x + Q+jO + 3jwxo)~ l (24) 

For the second-order approximation, we see that we get neither constant terms nor even terms in (±nj wt). 
This will remain true when iteration is carried out to higher order. We conclude that a solution based on 
iteration of the linear solution will be capable of representing 01 ly the first symmetric regime. 

Next, we shall discuss the validity of the harmonic balance method. The harmonic balance method 
is based on the truncated Fourier series: 

N N 

x(t) = Y^A n co3(nwt) + ^BnainCmvt) (25) 

n=0 n=0 

For the first-order approximation of a symmetric solution, we have: 

x,(t) = Aicoa(wt) + Bisin(wt) (26) 

When we substitute equation ( 26 ) into equation (1 1) and compute the coefficients A\ , B\ by matching the 
terms multiplied by cos{ wt) and sin( wt) and neglecting the higher-order terms multiplied by cos( mvt) 
and ain( mvt) with n > 1 , we obtain the solution via harmonic balance. Computation of the magnitudes of 
the neglected terms involves some algebraic manipulation that becomes unwieldy for higher-order approx- 
imations. The Laplace-Botel transformation has the merit of involving relatively easy algebraic manipu- 
lation that bears a strong analogy to the Heaviside operational calculus. Let us illustrate for the first-order 
approximation 

Xg > ' f - G(Xg) ( 27 ) 

Gi.x a ) = A- i(l — jwxo)~ l + A + i ( 1 + jwxq)~ x ( 28 ) 

where: 

Ati = ±jB x )/2 (29) 

If furthermore we neglect the transient terms, we have: 


G(x a ) = — 0(1 — o,x 0 )- 1 a:o(l - a 2 xo)- l G(x s )VG(x t )UG(xs) 

+ (1 — oia:o)“ 1 a:o( 1 — a2Xo)~ 1 xo(f/2)((l — jwxo)~ x + (1 + jwx o) -1 ) (30) 

Decomposition into simple rational fractions yields: 
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jivxo) 1 + An(l + jwxo)" 1 


= ^.{^ (U -ytyxo)' 1 + (l + /ti»a 0 )“ l ) -/?( 


»)-» 

(tt^— w 2 )+;'art« 


+ 


3^4_i ( 1 + jwxp ) + i4j|(l — .^tuxo)" 1 -4+i(l + B/toxo) -1 

(u>o — tu 2 ) + ;o(ty (wq — 9w 2 ) + /a3tu (tt>o — 9iu 2 ) — /a 3to 


( 31 ) 


Here the important point is that we can identify the amplitude ratio of the component of [lw] to that of 
[3w] through the factor /? 

- s I ( 32) 

(wo — u> 2 ) 2 + (ait;) 2 * (wo — 9tt> 2 ) 2 + (a3ty) 2 

with algebraic manipulation similar to that of the Heaviside operational calculus. Hence, we are able to 
find the parameter range in which the component in [3w] can be neglected. 


SYMMETRY-BREAKING BIFURCATION OF A HARD SPRING 


The idea underlying the study of symmetry-breaking bifurcation originates from the problem of 
predicting critical values of parameters at onset of chaos. The Melnikov criterion predicts onset of chaos 
in the case of a soft spring with some limited success (ref. 8); however it cannot be applied to the case of a 
hard spring. A natural alternative approach is then to predict critical values of parameters for the first few 
bifurcations, on the hypothesis that this may be sufficient to establish the trend toward the onset of chaos. 

It has been found (ref. 4) that the first-order approximation of the harmonic balance method cannot 
allow a symmetry-breaking bifurcation, while the second-order approximation may do so. Therefore, the 
validity of representing a solution by a Fourier series that is truncated at some order becomes questionable 
(ref. 4). There is no general theory that would allow us to affirm the validity of a truncated series solution. 
However, we shall show that a second-order harmonic balance representation leads to values of critical 
parameters at a symmetry-breaking bifurcation that are in reasonable agreement with numerical results 
based on a complete integration of Duffing’s equation. 

In order to visualize the symmetry-breaking bifurcation, we define the amplitude as < x 2 >- 
< x > 2 , where <> denotes the average over a period. We draw the amplitude response curve versus 
frequency of the symmetric solution of Duffing’s equation represented by 

(i) a first-order approximation 

x* !) = A\coawt + Bisinwt (33) 


(ii) the response curve of the symmetric solution represented by a second-order approximation 

x< 2) = Aj coswt + B\ simvt + A 3 cos 3 wt + B 3 stn3 wt ( 34) 
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and (iii) the response curve of the asymmetric solution represented by: 

Xas = A"q + A'I coswt + B"i simvt + A '2 cos 2 wt + B 2 siv2 wt + A^ cos3 wt + £3 sin3 wt (35) 
The principal steps in using the Laplace-Borel transformation are as follows. 

1. Calculate the shuffle product. 

2. Decompose into simple rational fractions. 

3. Match coefficients of terms ( 1 — j nwxo ) _1 , where n is an integer. 

4. Solve numerically for the coefficients A n ,B n . 

The results are given in figure 1 for the following values of parameters: ct = 0 .2 , tu 0 = 1 .0 , 

/= 10 . 0 . 

The response curve of the symmetric solution represented by equation (32) consists of only one 
main resonance peak, while the one represented by equation (33) consists of a main resonance peak R\ 
and a secondary resonance peak R $ . The response curve of the asymmetric solution represented by equa- 
tion (34) consists of only a resonance peak R 2 , at a frequency range between R\ and £3 • In contrast to the 
case of a soft spring where the asymmetric solution is present for a wide range of frequency w and even 
for small values of the forcing amplitude /, the asymmetric solution for the hard spring is confined to a 
small range of w. Furthermore, it does not appear for small values of /. The difference can be viewed as 
originating from the difference in form of the potential well for the hard spring and for the soft spring. In 
the former case, the potential consists of a single symmetric well with infinite amplitude at infinite distance, 
while in the latter case, the potential well is confined to finite values. For the infinite symmetric well, the 
frequency has to be tuned to some frequency in order to break the symmetry of the motion, even at sufficient 
values of /. A complete numerical computation would include more resonance peaks; however, the main 
features at frequency w > 0.6 of figure 1 are in qualitative agreement with numerical results of reference 6 . 
Using the second-order harmonic balance approximation, we study its symmetry-breaking bifurcation by 
stability analysis based on Floquet’s theory. Our approach incorporates the mathematical works on Hill’s 
equation (ref. 7). The results are given for comparison with numerical results from reference 4 in figure 2 
(Values of parameters: ac = 0.2, ft = (2) 3 / 2 /3). 

We note first that we are able to reproduce one band of values of parameter f (usual approaches 
[4,5] reproduce only a lower limit value for /). Our limitation to one band of values originates from the 
second-order approximation for the symmetric solution. The values of / obtained for low values of w 
iw < 0.9) are incorrect. This is expected because an exact solution would include more harmonics and 
the amplitude of oscillation at low w is then affected. The results obtained for / at high values of w are in 
reasonable agreement with numerical results. 
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CONCLUSIONS 


We have applied the Laplace-Borel transformation to discuss representation of analytical solu- 
tions of Duffing’s equation. Exploiting the symmetry property of this equation, we have demonstrated the 
following. 

1. In an iterative procedure, use of the first order iteration based on the linear equation leads only 
to a description of the first regime (symmetric solution). 

2. Implementing the harmonic balance method with the Laplace-Borel transformation gives an 
algebra very similar to that of the Heaviside operational calculus. 

Analysis of the nonlinear mechanics of a hard spring by means of the Laplace-Borel transformation 
gives values of parameters at a symmetry-breaking bifurcation in reasonable agreement with numerical re- 
sults. It should be possible to implement the Laplace-Borel transformation on computers as a combination 
of symbolic and numerical computation in order to predict values of parameters for the first few period- 
doubling bifurcations. 
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Figure 1.- Amplitude response curve versus frequency for Duffing’s equation: x+0.2x + x + x 3 

fcos(wt). 
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Figure 2.- Values of forcing amplitude / versus frequency w for bifurcation to symmetry-breaking. 
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